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Abstract 

An ideal Maxwell-Boltzmann gas confined in various rectangular nano domains 
is considered under quantum size effects. Thermodynamic quantities are calculated 
from their relations with partition function which consists of triple infinite sum¬ 
mations over momentum states in each direction. To get analytical expressions, 
summations are converted to integrals for macro systems by continuum approxi¬ 
mation which fails at nanoscales. To avoid both from the numerical calculation of 
summations and the failure of their integral approximations at nanoscale, a method 
which gives an analytical expression for single particle partition function (SPPF) is 
proposed. It’s shown that dimensional transition in momentum space occurs at cer¬ 
tain magnitude of confinement. Therefore, to represent SPPF by lower-dimensional 
analytical expressions becomes possible rather than numerical calculation of sum¬ 
mations. Considering rectangular domains with different aspect ratios, comparison 
of the results of derived expressions with those of summation forms of SPPF is 
done. Its shown that analytical expressions for SPPF give very precise results with 
maximum relative errors of around 1%, 2% and 3% at just the transition point for 
single, double and triple transitions respectively. Based on dimensional transitions, 
expressions for free energy, entropy, internal energy, chemical potential, heat ca¬ 
pacity and pressure are given analytically valid for any scale. 

Keywords: Quantum size effects; Dimensional transitions; Nano thermodynamics; 
Dimensional superpositions 


1 Introduction 

Advancements in nanoscience and nanotechnologies increased the importance of studies 
related to the nanoscale properties of matter, in recent years. On the presence of quantum 
conhnement effects, thermodynamic properties of systems differ from those of macro 
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scale. Due to wave character of particles, quantum size effects (QSE) appear at nanoscale 
systems as a consequence of quantum confinement. 

Thermodynamics under QSE became an active research area and many novel aspects 
such as quantum boundary layer, anisotropic gas pressure, loss of additivity in extensive 
properties, quantum forces and gas diffusion due to size difference in both Maxwell- 
Boltzmann (MB) and quantum (Fermi-Dirac and Bose-Einstein) gases have been studied 
in literature [hfjlT] . Even though size and shape dependencies in thermodynamics of nano 
systems has been recently studied, dimensional transitions due to QSE in thermodynamic 
quantities and transition conditions have not been examined in literature. 

Thermodynamic state functions are represented by summations from their dehnitions 
in statistical mechanics 18 ,19 . In thermodynamics of macro systems, it is common 


to replace summations by integrals to decrease the calculation burden and to obtain 
analytical expressions. Conversion of summations to integrals is done by using continuum 
approximation and density of states concept, under the assumption of discrete energy 
levels are so close to each other that they can be considered as continuous variables. 
Continuum approximation works extremely well in macro scale. However, in nanoscale, 
ascending divergence appears between summations and their integral approximations 

00 


At meso scales, instead of integral approximations, Weyl’s conjecture or hrst two 
terms of Poisson summation formula (PSF) are used in literature §1131 [^[^. On the 
other hand, when characteristic sizes of the domain are smaller than half of the most 
probable de Broglie wavelengths of particles, even these approximations considered in 
literature become insufficient. In the case of such strong conhnements, one should either 
consider the numerical calculation of exact summations in thermodynamic properties or 
look for a better solution. 

In this article, a method based on dimensional transitions in momentum space due to 
quantum conhnement is proposed for the calculation of single particle partition function 
(SPPF) of MB statistics. It is seen that the method leads to analytical expressions which 
give nearly the same results of exact summations. Then from the analytical expressions 
of SPPF, thermodynamic quantities such as free energy, entropy, internal energy, heat 
capacity and pressure are derived analytically. Consequently, the results successfully 
represent thermodynamic behavior of strongly conhned systems without need to calculate 
relevant triple summations. 


2 Determination of Dimensional Transition Point in 
Momentum Space 

For a high temperature and/or low density ideal gas, SPPF is given in triple summation 
form by using MB statistics as follows 

c = E-p(-(^) w 

where in = {1,2,3,...} are momentum state variables with subscript of n = {1,2,3} 
indicating directions, £ are the energy eigenvalues from the solution of the Schrodinger 
equation for a rectangular domain, ks is Boltzmann’s constant and T is temperature. By 
dehning a conhnement parameter = L^Ln that indicates the magnitude of conhne- 
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ment in n’th direction, Eq. (1) can also be written as follows 


3 oo 

C = n 

n=l 2n = l 


( 2 ) 


In the definition of a, Ln is the size of the domain in the n’th direction and Lc = 
h/y/SmkBT is the half of the most probable de Broglie wavelength of particles where h 
is Planck’s constant and m is mass of the particle. 

Although triple summations are the exact way to calculate SPPF, in macro scales 
where a values are much smaller than unity {a << 1), summations are replaced by 
integrals making continuum approximation. On the other hand, when domain size in 
particular direction (L„) is on the order of Lc, like in nanoscale {a ~ 1), there is a 
considerable deviation between the results of summations and integrals. Therefore, a 
more precise evaluation method, such as Poisson summation formula (PSF), is necessary 
to calculate summations. PSF is written for even functions, f{i) = f{—i), as 


i=i do 
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By evaluating each term of PSF, single particle partition function in ID can be exactly 
represented as 


CiD = 5]] exp 


*1=1 
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(4b) 


The hrst term of PSF is the conventional integral term and it’s accuracy for macro 
domains (a —)■ 0) is quite well. When domain size gets closer to Lc (a —)■ 1), contribution 
of the second term (called zero correction term) of PSF becomes appreciable. First two 
terms of PSF can be called one-dimensional (ID) representation, since it represents the 
contributions of whole momentum state variables, A = {!,..., cx)}. Contribution of the 

MEHH 


zero correction term is studied as QSE in literature 1708 


especially for weakly 


conhned domains (a < 1) in which the contribution of the third term of PSF is negligible. 
To fully represent SPPF in any scale, even for strongly conhned domains (a > 1), whole 
terms of PSF have to be considered. However, since there is no known analytical solution 
for the third term which include another inhnite summation, either we have to calculate 
inhnite summations of SPPF numerically or End a method to give analytical expressions 
representing the exact behavior of SPPF successfully. For strongly conhned cases («>!), 
fortunately, contribution of just the hrst term (A = 1) of SPPF’s summation becomes 
extremely dominant. Therefore, it can be called as zero-dimensional (OD) representation 
since it represents almost only the contribution of the ground state. 

To see the success of ID and OD representations, variations of their analytical expres¬ 
sions with conhnement parameter are compared with that of exact ID SPPF based on 
numerical calculation of summation, in Figure 1. The solid curve shows the variation 
of exact form of ID SPPF (Eq. 4a) with conhnement parameter ai. Plus signs show 
the results of ID analytical representation based on hrst two terms of PSF, while cross 
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signs demonstrate the results of OD representation based on just the hrst term = 1) 
of SPPF’s summation. 



ai 


Figure 1: Dimensional transition from ID to OD representation in SPPF for an MB gas 


It is understood that ID analytical expression matches perfectly with the exact form 
of SPPF until the certain value of a = a*, after that it considerably deviates. In the same 
manner, OD expression substantially deviates before a*, while it perfectly matches with 
the solid curve after that. This certain value of a, which corresponds to the transition 
point from ID to OD representation, can be found as a* = 1.23547... by equating ID 
expression to OD one. Therefore, an analytical expression representing SPPF for any 
values of a can be given as follows 


OO y , -X C 

^exp f [exp(-a2)]^-" 

il=l 


c J c = 1 for ai < a* 
2 ] c = 0 for ai > a* 


(5) 


where the transition value is a* = 1.23547.... This approximation gives 1% error even 
at the transition point. The transition represents the dimensional change (contrac¬ 
tion/expansion) of momentum space in particular direction. Henceforward, the domain 
in the n’th direction will be considered as free and conhned for < a* and > a* 
respectively. 


2.1 Dimensional Transitions and Superpositions in Single Par¬ 
ticle Partition Function 

By using the dimensional transition point in momentum space, we may now consider 
SPPFs for the domains which are conhned at nanoscale in certain directions. 3D domain 
is free in all three directions, so that conhnement values in all directions are smaller than 
a*, the transition value. For 2D domain, ai and for ID domain, ai and 02 are chosen as 
greater than a* respectively. OD domain is considered as conhned in all directions, which 
means all of their conhnement parameters are larger than a*. In this sense, using Eq. 
(5), we can write analytical expressions for SPPF of an ideal MB gas for 3D, 2D, ID and 
OD domains as 
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^CoD = exp [-al - al - al) , {ai, 02 , 03 } > a* 


Although dimensional transition point has a unique value, by engaging the sizes of 
the domain to each other in each direction, it is possible to obtain different transition 
points correspond to different dimensional transitions for anisometric domains. In order 
to engage sizes to each other, we may define constant aspect ratios between sizes in three 
directions. Therefore, we can change the size of the domain by changing only one of the 
sizes, which is chosen here as Li, while preserving its shape. Aspect ratio between sizes 
in first and second directions is denoted by ri 2 = T 1 /L 2 and while the second and third 
directions by r 23 = L 2 /T 3 . In that case, confinement parameters of second and third 
directions can be written in terms of the confinement parameter in the first direction as 
0^2 = Q!iri 2 and = airi 2 r 23 respectively. In Figure 2, relative errors of analytical 
expressions given by Eq. ( 6 ) are given for different aspect ratios. 






Figure 2: Relative errors of analytical expressions of SPPF for four various aspect ratio 
sets. Peaks correspond to the relevant dimensional transition points, a) Fully anisometric 
domain with aspect ratios ri 2 = r 23 = 0.7. b) Semi-anisometric domain with aspect 
ratios ri 2 = 1 and r 23 = 0.7. c) Semi-anisometric domain with aspect ratios ri 2 = 0.7 
and r 23 = 1. d) Isometric domain with aspect ratios ri 2 = r 23 = 1. 

Peaks in Figure 2 correspond to the values a*, a*/r 12 and a*/ri 2 r 2 ^ for transitions from 
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3D-2D, 2D-1D and ID-OD respectively. As it is seen from Figure 2a, each dimensional 
transition causes an error of 1% even at the transition points. If the sizes in two different 
directions are the same, dimensional transitions in these directions superpose and it 
causes 2% error as it is seen in Figures 2b and 2c. Similarly, for an isometric domain, 
superposition of three different dimensional transitions cause 3% error. Figure 2d. 


3 Dimensional Transitions in Thermodynamic Prop¬ 
erties 


Free energy of an ic 

6 cl 

V 

B gas consisting of N particles is simply given by P = 

. Dimensionless free energy can be dehned as P = F/NksT. 

-Nk,T[\n{C/N) +1] 

18 
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Then, we can obtain analy 

deal expressions for dimensionless free energies of a MB gas 


for 3D, 2D, ID and OD cases by considering dimensional transitions 


Fsd = -In 
F2d = -In 
Fid = -In 
Fqd = -In 

respectively, where qn = «n/ — ^n) with n = {1,2,3} as direction index. It is clear 

that ttn is always smaller than a/F since Qn is necessary only for the case of < a*. 
Note that, we can use MB statistics if and only if 1 << N « l/aia 20 i^ condition holds. 
This condition holds even if one of the conhnement parameters an is smaller enough than 
unity. Because in that condition the system contains enough particles to do statistics 
{N >> 1) and interparticular distance is much greater than to use MB statistics. 
Since OD case does not satisfy this condition, one could not use MB statistics for OD 
conditions. Nevertheless, in this article, we gave also OD expressions of thermodynamic 
properties, just to give a mathematical intuition. Therefore, it is worth to note that 
thermodynamic expressions given for OD do not represent a meaningful physical situation 
since the condition of N « l/aia 2 a^ leads N < 1 for OD condition. 

From free energies, we can derive all other thermodynamic quantities such as entropy, 
internal energy, chemical potential, heat capacity at constant sizes and pressure as 


maia 2 aj “ ^ + 1 ) + (^2 + 1) + In + 1 ) 


TT 


ANa2as J 
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2Na-^. 


— 1 + + /n (q '2 + 1) + /n {q^ + 1), { 0 ^ 1 } > (F 

— 1 + a^ + a\ + In [q^ + 1 ), {q;i, ^ 2 } > a* 
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— ) — 1 + + ^2 + «3, {ai,a2, as} > a* 
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/ N,V 

C=(^-^ 
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U = F + TS 
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Li,L2,L3 


P =^kT^ 

"" C \dLn 


Thermodynamic properties can be written in their dimensionless forms as S = S/Nks, 
U = U/NkeT, J1 = n/NkeT, C = C/Nke and An = PunV/NkeT where V = L 1 L 2 L 3 . 
Even though pressure is a scalar quantity in macro scale, it becomes a tensorial quantity 
with zero off-diagonal elements for anisometric domains. In other words, Pn 7 ^ P22 ^ P33 
for an anisometric rectangular domain where Li ^ L 2 ^ T 3 . 
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Dimensionless entropy, internal energy, chemical potential, heat capacity and pressure 
expressions for MB gases are derived by considering dimensional transitions and given in 
Table 1. 

Table 1; Thermodynamic properties of Maxwell-Boltzmann gases in various dimensional 

conditions _ 

Entropy 

S3d = 1 + In ~ (^1 + l)-ln (gs + 1) - In (gs + 1) + + ^2 + gs) 

S 2 D = 2 + In 4 ^^ 203 ) ~ ^^2 + 1) — In (gs + 1 ) + | (g 2 + gs), {«i} > 

SiD = l + In ( 2 ^) - (?3 + 1 ) + Kgs), {tti, 02 } > a* 

Sod = ^ + In , {ai, as, as} > a* 

Internal energy 

U 3 D = I + ^ (gi + g2 + gs) 

h^ 2 D = 1 + «! + ^ (g 2 + gs) ) {^ 1 } > 

UiD = I + af + af + |(gs), {cti, as} > a* 

Uqd = a^ + a^ + a|, {cii, cis, as} > a* 

Chemical potential 

/fsD = -In (sTv^^^) + /n (gi + 1 ) + /n (g 2 + 1 ) + /n (gs + 1 ) 
h2D = —In ^47Va2 03 ) + (g2 + 1) + /u (gs + 1), {ai} > a* 

JliD = -In (^2^) +al + al + ln (gs + 1), {ai, a 2 } > a* 

JioD = -In + al + al + a|, (ai, a 2 , as} > a* 

Heat capacity at constant sizes {Li,L 2 ,Ts} 

Cod = 1 + 1 [gi(i “ gi) + g2(i ~ 0.2) + gs(i ~ gs)] 

C 2 D = 1 + i [g 2 (i - g 2 ) + gs(i - gs)], {ai} > a* 

CiD = 1 + 1 [gs(i - gs)], {tti, ^ 2 } > a* 

Cod = 0, {ai,a 2 ,as} > a* 

Pressure 

p3/^ = (gi + l), = + + 

Ph^ = 2al Pif = {q 2 + 1), Pii^ = {qo + l), {a,} > a* 

Pn = 2tt?, ^ 22 ^ = 2ai, P^f = (gs + 1), {ai, a 2 } > a* 

P^r = 2al Pl^ = 2al Pl^ = 2al {ai, a2, as} > a* 
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It is easy to check that, all expressions derived in this article go to their classical 
expressions when conhnement parameter q;„’s go to zero. So they are valid not only for 
nano systems but also for meso and macro systems. They are generic in that sense, within 
the error limit at dimensional transition points. 

Note for heat capacity that, instead of referring to constant volume, we refer here 
constant sizes which means size of each direction Li, L 2 and L 3 are constant. Otherwise, 
combination of different Ln values can also correspond to the same volume, although they 
correspond different heat capacity values, different from that in macro scale. 

Unlike the scalar behavior of pressure in macroscopic scale, quantum confinement in 
one direction brings a positive contribution to the pressure in that direction which make 
the pressure a tensorial quantity in nanoscale. Besides, when conhnement of a particular 
direction is larger than a*, pressure at that direction grows quadratically. 


4 Conclusion 

By introducing dimensional transitions of thermodynamic properties for MB gases, it’s 
shown that it is possible to obtain analytical expressions within a negligible error to rep¬ 
resent the exact forms of thermodynamic properties based on triple inhnite summations. 
In this regard, analytical results are obtained for SPPF, free energy, entropy, internal en¬ 
ergy, chemical potential, heat capacity and pressure, for 3D, 2D, ID and OD cases. It is 
seen that thermodynamic quantities have strict size dependency which cause dimensional 
contraction in free momentum space for nano systems. 

Critical value of conhnement parameter is found as a* = 1.23547.... When the conhne¬ 
ment parameter in a particular direction becomes greater than critical value of transition, 
momentum component in that direction can only take its ground-state value. Then, de¬ 
gree of freedom of momentum space decreases one unit. Therefore, depending on the 
values of conhnement parameter in each direction, degree of freedom takes the values of 
3, 2,1 and 0. By considering this dimensional transition model, it is possible to express 
the thermodynamic properties analytically with 1 % error for any conhned system. 

As it is expected, classical expressions can be recovered in case of the absence of 
conhnement, (q;„ —)■ 0). Dimensional transitions in ideal Fermi and Bose gases are under 
consideration. Thermodynamic properties of strongly conhned ideal Fermi and Bose gases 
in nano domains will then be analytically obtained for any scale. 
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